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ABSTRACT 

When Earth-mass extrasolar planets first become detectable, one challenge will be to determine 
which of these worlds harbor liquid water, a widely used criterion for habitability. Some of the first 
observations of these planets will consist of disc-averaged, time-resolved broadband photometry. To 
simulate such data, the Deep Impact spacecraft obtained light curves of Earth at seven wavebands 
spanning 300-1000 nm as part of the EPOXI mission of opportunity. In this paper we analyze disc- 
integrated light curves, treating Earth as if it were an exoplanet, to determine if we can detect the 
presence of oceans and continents. We present two observations each spanning one day, taken at 
gibbous phases of 57° and 77°, respectively. As expected, the time-averaged spectrum of Earth is 
blue at short wavelengths due to Rayleigh scattering, and gray redward of 600 nm due to reflective 
clouds. The rotation of the planet leads to diurnal albedo variations of 15-30%, with the largest 
relative changes occuring at the reddest wavelengths. To characterize these variations in an unbiased 
manner we carry out a principal component analysis of the multi-band light curves; this analysis 
reveals that 98% of the diurnal color changes of Earth are due to only 2 dominant eigencolors. We use 
the time- variations of these two eigencolors to construct longitudinal maps of the Earth, treating it as 
a non-uniform Lambert sphere. We find that the spectral and spatial distributions of the eigencolors 
correspond to cloud-free continents and oceans; this despite the fact that our observations were taken 
on days with typical cloud cover. We also find that the near-infrared wavebands are particularly useful 
in distinguishing between land and water. Based on this experiment we conclude that it should be 
possible to infer the existence of water oceans on exoplanets with time-resolved broadband observations 
taken by a large space-based coronagraphic telescope. 

Subject headings: methods: data analysis — (stars:) planetary systems — 



1. INTRODUCTION 

The rate of discovery of extrasolar planets is increas- 
Lg and every year it is possible to detect smaller plan- 
;s. It is only a matter of time before we detect Earth- 
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analogs, but even then our ability to study them will 
remain limited. Due to exoplanets' great distance from 
us and their relative faintness, spatially resolving them 
from their host stars is only currently poss ible for hot, 
youn g Jovian plane t s in l o ng-period orbits (|Kalas et al.l 
120081 iMarois et all l2008t iLagrange etlHI l2008f ). Fot 
other planets — including Solar System analogs — spa- 
tially separating the image of the planet from that of the 
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host star will have to wait for space-based telescopes like 
TPF/Da rwin ()Traub et alJ 120061 : iBeichman eFall l200d 
lFridlundl[200l . But even these telescopes will not have 
sufficient angular resolution to spatially resolve the disc 
of an exoplanet. 

As noted over a century ago bv iRusselH (|1906( l. vari- 
ations in the reflected light of an unresolved rotat- 
ing object can be used to learn about albedo mark- 
ings on the body. Such light curve inversions have 
proved valuable to interpret the photometry of ob- 
jects viewed near full phase and led, for example , 
to the first albedo map of Pluto (|Lacis fc Fixl HqtI ). 
More recently, thermal light curves have made it 
possible to measure the day/n ight temperature con- 
trast of sh ort-period cxoplancts (Harrington et al1l2006l : 
iCowan et al. 2007 : Snellen et al. 2009). Light curve in- 
version (jCowan fc A gol 2008|) has even been used to con- 
struct coarse longitudinal t hermal maps of hot Jupiters 
(|Knutson et al . I [20071 . [2009? ) . 

The optical and near-IR light curves of Earth, on the 
other hand, have not been thoroughly studied to date. 
Earthshine, the faint illumination of the dark side of 
the Moon due to reflected light from Earth, has been 
used to study the reflectance spectrum, cloud cover vari- 
ability, vegetation signatures and the effects of specular 
reflection for limited region s of our planet ([Goode et ahl 
2001; Woolf ct al. 2002; Qiu et al.ll2qollPalle et al.ll200l 
2004 Montan es-Rodriguez'et al .' '20051: ISeager et al 
2005[ Hamda ni et all 120061: iMontancs-Rodriguez et al 



20061 I2007t iLangford et all [20b9i ). Brief snapshots 



of Earth obtained with the Galil eo spacecraft have 
been used to stud y our planet (jSagan et al.l 119931 : 
iGeissler et al.l Il995| ) and numerical models have been 
developed to predict how diurnal variations in disc- 
integrated lig ht c o uld be used t o char act erize Earth 
(Ebrd ct al. 200l iTinetti et a p l2006airbl : iPalle et all 
[20081: IWilfiams fc Gaidojl2008h . 

This paper is the first in a series analyzing the photom- 
etry and spectroscopy of Earth obtained as part of the 
EPOXI mission and is written in a different spirit than 
most studies of Earthshine. Rather than attempting to 
produce a detailed model which exactly fits the observa- 
tions (a.k.a. "forward modeling" ), we make a few reason- 
able simplifying assumptions which allow us to extract 
information directly from the data ("backward model- 
ing" ) . This approach is complementary to detailed mod- 
eling and will be especially appropriate when studying an 
alien world with limited data. This paper is organized 
as follows: in § 2 we describe the time-resolved obser- 
vations of the entire disc of Earth used in this study; in 
§ 3 we use principal component analysis to determine the 
dominant spectral components of the planet in a model- 
independent way; we use light curve inversion in § 4 to 
convert the diurnal albedo variations into a longitudinal 
map of Earth; we discuss our results in § 5; our conclu- 
sions are in § 6. 

2. OBSERVATIONS 

The EPOXI"' mission reuses the still-functioning Deep 
Impact spacecraft that successfully observed comet 

The University of Maryland leads the overall EPOXI mission, 
including the flyby of comet Hartley 2. NASA Goddard leads the 
exoplanet and Earth observations. 



9P/Tempel 1. EPOXI science targets include several 
transiting cxoplancts and Earth en route to a flyby of 
comet 103P/Hartley 2. The EPOXI Earth observations 
are valuable for exoplanet studies because they are the 
first time-resolved, multi-waveband observations of the 
full disc of Earth. These data reveal Earth as it would 
appear to observers on an extrasolar planet, and can only 
be obtained from a relatively distant vantage point, not 
from low-Earth orbit. The data consist of an equinox 
observation on March 18, and near solstice on June 6. 
An observation taken on May 28 included a planned lu- 
nar transit, but we save the analysis of those data for a 
future paper. The observations, summarized in Table [H 
were taken when Earth was in a partially illuminated — 
gibbous — phase. With the rare exception of transiting 
planets, this is the phase at which habitable cxoplancts 
will be observed. 

Deep Impact's 30 cm diameter t elescope coupled with 
the High Resolution Imager (HRI, [Hampton et al.l[2005| ) 
recorded images of Earth in seven 100 nm wide opti- 
cal wavebands spanning 300-1000 nm. Hourly observa- 
tions were taken with the filters centered on 350, 750 and 
950 nm, whereas the 450, 550, 650 and 850 nm data were 
taken every 15 minutes; each set of observations lasted 
24 hrs. The exposure times for the different wavebands 
are: 73.4 msec at 350 nm; 13.3 msec at 450 nm; 8.5 msec 
at 550 nm; 9.5 msec at 650 nm; 13.5 msec at 750 nm; 
26.5 msec at 850 nm; 61.5 msec at 950 nm. Although the 
EPOXI images of Earth offer spatial resolution of better 
than 100 km, we mimic the data that will eventually be 
available for cxoplancts by integrating the flux over the 
entire disc of Earth and using only the hourly EPOXI ob- 
servations from each of the wavebands, producing seven 
light curves for each of the two observing campaigns, 
shown in Figure [T] Our results are the same when we 
use the 450, 550, 650 and 850 nm data from :00, :15, 
:30 or :45. The photometric uncertainty in these data is 
exceedingly small: on the order of 0.1% relative errors. 

Since we are interested in the properties of the planet 
rather than its host star, we normalize the light curves 
by the average solar flux in each bandpass using the so- 
lar spectrum^ to obtain the reflectivity in each wave- 
band. We express the brightness of the planet as an 
apparent albedo, the average albedo of regions on the 
planet that are both visible and illuminated during each 
observation^. The spacecraft was above the middle of 
the Pacific Ocean at the start of both observing cam- 
paigns, so the shapes of the light curves are similar for 
both epochs. The June light curves vary more rapidly 
because a smaller fraction of the illuminated hemisphere 
of Earth was visible (62% rather than 77%) and thus less 
of the planet was averaged together in any given frame. 

2.1. Cloud Variability 

Clouds coyer ro ughly half of Earth at any point in time 
([Palle et al.[[2008D and th ey dominate the disc -integrated 
albedo of the planet (e.g.. [Tinetti et al][2006a| ). Mapping 

5 We use the ASTM-E-490 Standard Solar Constant 
and Z ero Air Mass Solar Spectral Irradiance Tables: 



|http: / / www. astm.org/Standards/E490. htm 

° The ratio of the observed flux to the expected flux for a planet 
of the same size and phase exhibiting diffuse reflection and with an 
albedo of unity every where on its surface. Mathematical details of 
this definition are in i; l4.1l 
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TABLE 1 

EPOXI Earth observing campaigns 



Date 


Starting 


Phase 


Illuminated Fraction 


Spacecraft 


Angular Diam. 


Pixels Spanned 




CML 




of Earth Disc 


Range 


of Earth 


by Earth 


3/18/2008 


150° W 


57° 


77% 


0.18 AU 


1.63' 


~ 240 


5/28/2008 


195° W 


75° 


63% 


0.33 AU 


0.89' 


130* 


6/4/2008 


150° W 


77° 


62% 


0.34 AU 


0.87' 


~ 130 



The CML is the Central Meridian Longitude, the longitude of the sub-observer point. The 
planetary phase, a, is the star-planet-observer angle and is related to the illuminated fraction by 
/ = ^(1 + cos a). *The 5/28/2008 observation is not used in this paper due to a planned lunar 
transit. 
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Fig. 1. — Seven light curves obtained by the EPOXI space- 
craft on March 18 (solid lines) and on June 4, 2008 (dashed lines). 
The bottom-right panel shows changes in the bolometric albedo of 
Earth. 

surface features could be problematic if large scale cloud 
formations move or disperse on timescales shorter than 
a planetary rotation. 

The necessary condition for variable cloud patterns — 
surface pressure and temperature near the condensation 
point of water — is a likely precondition for habitability, 
and hence may pose a problem for the planets that in- 
terest us the most. Changeable cloud cover may indicate 
the presence of water vapor in a planet's atmosphere, but 
here we are concerned with the presence of liquid water 
on the planet's surface. 

After 24 hours of rotation the same hemisphere of 
Earth should be facing the Deep Impact Spacecraft so 
the integrated brightness of the planet's surface should 
be identical, provided one has accounted for the dif- 
ference between the sidereal and solar day, as well as 
slight changes in the geocentric distance of the space- 
craft and in the phase of the planet as seen from the 
spacecraft. Note that our observations are not taken 
near f ull phase so we ma y safely ignore the opposition 
effect (|Hapke et al.l 119931 ) . Even after correcting for all 
known geometric effects, the observed fluxes at the start 
and end of a given observing campaign differ by 2.2% 
and 3.4% for the March and June observing campaigns, 
respectively, as shown in Figure [2] We attribute this dis- 
crepancy to diurnal changes in cloud cover. Therefore, 
even though the EPOXI photometry is excellent, for the 
purposes of our model fits we use effective uncertainties 
equal to |Fstart — ^Gnd|/2 for each waveband. 

Our 24-hour cloud variability of a few percent is some- 
what small er than estimates from Earthshine o bserva- 
tions (e.g.. iGoode et al.l l200lt iPalle et all |2004 found 



Fig. 2. — The discrepancy in apparent albedo between the start 
and end of each observing campaign. Note that during the March 
observations the flux decreased at all seven wavebands, while in 
June the flux increased. We attribute these changes in albedo to 
changes in cloud cover. 

day-to-day cloud variations of roughly 5% and 10%, re- 
spectively), and is a small effect compared to the rota- 
tional modulation of Earth's albedo or the 10-20% dif- 
ferences between the light curves from the two EPOXI 
observing campaigns. The 10-30% variations in apparent 
albedo we observe due to the Earth's rotation (with the 
largest variations occuring at near inf rared wavebands) 
agrees with previ ous optical studies (jFord et al.l 120011 : 
iGoode et al.ll200ll . find diurnal variations of 15-20%). 

The differences between the March and June obser- 
vations could be due to some combination of stochas- 
tic changes in cloud cover, coherent (seasonal) changes 
in cloud cover, or simply a change in viewing geometry 
(which we discount in § 4.3). Although daily changes in 
cloud cover are modest, the cloud cover will be entirely 
different for observa tions taken mont hs apart (see, for 
example, Figure 6 of iPalle et al]|2008f ). 

3. DETERMINING PRINCIPAL COLORS 

In this study we assume no prior knowledge of the dif- 
ferent surface types of the unresolved planet. Our data 
consist of 50 broadband spectra of Earth (25 hourly ob- 
servations for each of two epochs). The time-averaged 
spectrum of Earth is blue at short wavelengths due to 
Rayleigh scattering and gray longward of 550 nm be- 
cause of clouds, as shown in Figure [H The changes in 
color of Earth during our observations can be thought of 
as occupying a 7-dimensional parameter space (one for 
eac h waveband). P rincipal component analysis (PCA, 
e.g. IConnollv et"al|[l995 l allows us to reduce the dimen- 
sionality of these data by defining orthogonal eigenvec- 
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Fig. 3. — Time-averaged broadband spectrum of Earth, based 
on EPOXI observations taken on 03/18/2008 and 06/04/2008 (the 
time-averaged spectra for the two epochs are indistinguishable). 
The error-bars show la time-variability. The steep ramp at short 
wavelengths is due to Rayleigh scattering. The near-IR wavebands 
exhibit the largest relative time-variability. 

tors in the parameter space (eigencolors). Qualitatively, 
the observed spectrum of Earth at some time t can be 
recovered using the equation: 

7 

A*{t) = {A*)+Y,C,{t)A, (1) 

4 = 1 

where (A*) is the time-averaged spectrum of Earth, Ai 
are the seven orthogonal eigencolors, and d are the in- 
stantaneous projections of Earth's colors on the eigen- 
colors. The terms in the sum are ranked by the time- 
variance in Ci, from largest to smallest. Insofar as the 
projections are consistently small for the highest i, they 
can be ignored with only a minor penalty in goodness- 
of-fit. 

We find that 98% of the changes in color do not oc- 
cupy the whole parameter space but instead lie on a two- 
dimensional "plane" defined by the two principal eigen- 
colors. That is to say, truncating the sum in Equation [1] 
at i = 2 only leads to errors of a couple percent. In detail, 
the plane has some thickness to it: two additional com- 
ponents only present at the l-2a level that we neglect in 
the remainder of this study. We estimate the uncertainty 
in the PCA by creating 10,000 versions of the light curves 
with added Gaussian noise. The standard deviation in 
the resulting PCA parameters provides an estimate of 
their uncertainty. The primary eigencolors {Ai and A2 
from Equation[T]) are shown in Figure|3]and their relative 
importance as a function of time (Ci(i) and C2{t) from 
Equation [T]) is shown in Figure [S] The uncertainties on 
the PCA are correlated but we represent them as error 
bars in the figures. 

The eigencolors should not be thought of as spectra 
of different surface types. Rather, they are particular 
combinations of filters which are most sensitive to the 
different surface types rotating in and out of view. As 
such, the eigencolors are a relative color from the Earth 
mean. The first eigencolor is most sensitive to variations 
in the red wavebands (since the albedo of Earth varies the 
most at near-IR wavelengths) and the second eigencolor 
is most sensitive to the blue wavebands. 

Since the mean Earth spectrum is — to first order — 
a cloud spectrum seen through a scattering atmosphere. 
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Fig. 4. — Spectra for the two dominant eigencolors of Earth, 
as determined by PCA (top panel). For comparison, the bottom 
panel shows actual spectra of clouds, soil and oceans on Earth 
iTinetti et al.ll2006al: IMcLinden et al.|[T997l) . 
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Fig. 5. — Contributions of Earth's two principal colors, as deter- 
mined by PCA, relative to the average Earth spectrum from both 
epochs of observation. Each set of observations spans a full rota- 
tion of the planet, starting and ending with the spacecraft directly 
above 150° W longitude, the middle of the Pacific Ocean. 

the principal colors of Earth are related to the main sur- 
face types on Earth: cloud-fre e continents are mo st re- 
fiective at longer wavelengths (jTinetti et al.|[2006al con- 
structed from the ASTER Spectral Library found at 
http://speclib.jpl.nasa.gov), while cloud-free ocea ns are 
most reflective in the blue (McLinden et al.lll997| ). For 
example, the presence (or lack) of continents shows up 
as positive (or negative) excursions of the red eigencolor. 
The relative contributions of the colors are the projection 
of the Earth's instantaneous color onto an eigencolor. 
They are not identical in the March and June observa- 
tions due to different cloud cover: when both the red and 
blue eigencolors are positive, there was more than aver- 
age cloud cover, while regions with uniformly low eigen- 
colors correspond to relatively cloud- free regions. Never- 
theless, the similar shapes of the two sets of curves would 
indicate to extrasolar observers that the principal colors 
are most sensitive to fixed surface features that were vis- 
ible from one season to the next (We run the PCA on 
both the March and June observations simultaneously.). 
The implicit "third surface" in this analysis is the time- 
averaged spectrum of the Earth, which includes clouds 
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and Rayleigh scattering. 

To test how sensitively the results of the PCA depend 
on photometric uncertainty, we repeat the analysis with 
additional noise. Our principal result — the significant 
detection of red and blue eigencolors and their temporal 
variations — are essentially unchanged for photometric 
uncertainties smaller than 2-3%. Observations of exo- 
planets of this quality are not around the corner, but can 
be obtained in the foreseeable future: a 16 m space tele- 
scope (e.g., iPostman fc ATLAST Concept Studv Teaml 
|2009() equipped with a coronagraph could obtain 2% pho- 
tometry with 1 hour exposures of an Earth-analog at 
10 pc. We have assumed a planet with the same ra- 
dius and mean albedo as Earth, orbiting at 1 AU from 
a Sun-like star and obser ved ne ar quadrature. We com- 
pute signal-to-noise as in lAgol |2007) including photon 
counting noise from the planet and PSF noise from the 
host star, but neglecting zodiacal and exozodiacal noise. 

4. MAPPING SURFACE TYPES 

The time- variation in the eigencolors. Figure [5l tells 
us about the spatial variations of the colors around the 
planet. A region on the planet contributes more or less to 
the disc-integrated light depending on the amount of sun- 
light the region is receiving, its projected area as seen by 
the observer, and its albedo. As time passes, different re- 
gions of the planet rotate through the "sweet-spot" where 
the combination of illumination and visibility is optimal. 
The zeroth order approach to mapping the planet would 
be to assume that all of the light from the planet orig- 
inates from this region, as would be the case for purely 
specular reflection. The observed light curve could then 
be directly converted into a longitudinal albedo map of 
the planet at the appropriate latitude. We explore this 
limiting case in Appendix III. 

More exactly, determining the spatial distribution of 
a color based on the planet's multi-band light curves 
is a deconvolution problem equivalent to mapping the 
albedo markings on a body based on its disc-int egrated 
reflect ed light curve. This problem was solved bv lRussell 
(|1906l ) for outer Solar System objects, which are always 
observed near full phase as seen from Earth. Exoplanets, 
however, appear at a variety of phases (E.g., crescent, 
gibbous) and therefore require a more complete solution. 
We keep the problem tractable by considering the dif- 
fusely reflecting (Lambertian) regime, in which surfaces 
reflect light equally in all directions. 

Most materials are not perfectly Lambertian, instead 
scattering light preferentially backwards, forwards, or 
specularly. Coherent back-scattering is only signifi- 
cant near full phase (the opposition effect. iHapke et al 



19931) . In fact, o bservations of Earthshinc ( Qi u et al 
2OO3I lPalleetal.1 [200 3) and simulations (For d et al 



2OOH IWilhams fc GaidosI I200I . indicate that the disc- 
integrated light of a cloudy planet like Earth is well de- 
scribed by diffuse reflection, provided the star-planet- 
observer angle is close to 90°. Future missions — 
interferometers, coronagraphs, or occulters — will ideally 
observe planets at phase angles slightly smaller than 90° 
(a.k.a. "quadrat ure" ) , since that is when the S/N ratio 
is greatest (E.g.. I Agoll 120071 ). 

4.1. Reflected Light From A Non-Uniform Lambert 
Sphere 



The flux from a diffusely reflective non-uniform sphere 
can be parametrized in terms of its albedo map, A{9, 0), 
where 9 and are latitude and longitude on the planet, 
respectively. The visibility and illumination of a region 
on the planet at time t are denoted by V{d,(l>,t) and 
I{d,(f>,t), respectively. V is unity at the sub-observer 
point, drops as the cosine of the angle from the observer 
and is null on the far side of the planet from the observer; 
/ is unity at the sub-stellar point, drops as the cosine of 
the angle from the star and is null on the night-side of 
the planet. The mathematical details of V and / can 
be found in Appendix I. The planet/star flux ratio, e, is 
obtained by integrating the product of visibility, illumi- 
nation and albedo over the planet's surface: 

where Rp is the planet's radius and a is its mean orbital 
radius. The integral is over the entire surface of the 
planet, but the integrand is only non-zero for the regions 
of the planet that are both visible and illuminated (e.g., 
for -J of a planet viewed at quadrature). 

The flux ratio primarily depends on the planet's or- 
bital phase, the observer-planet-star angle, and the ra- 
tio (Rp/a)^. We define the apparent albedo. A*, as the 
ratio of the flux from the planet divided by the flux we 
would expect at the same phase f or a perfectly refl ecting 
{A = 1) Lambert sphere (see also lQiu et al.l[2003[ ): 



A*{t) 



VLAdn 

fvLdn ' 



(3) 



A uniform planet would have an apparent albedo that 
is constant over a planetary rotation; a true Lambert 
sphere would further have a constant apparent albedo 
during the entire orbit. For non-transiting exoplanets, 
the planetary radius is unknown, so A* can only be de- 
termined to within a factor of R^. 

4.2. Sinusoidal and N-Slice Maps 

The map A(0, 0) can take any form but diurnal 
brightness variations are due entirely to the longitude- 
dependence of albedo, provided that the planet's rota- 
tional period is much shorter than its orbital period. In 
the opposite extreme of a tidally-locked planet, seasonal 
variations in reflected light will be due to both albedo 
markings on the planet and changes in phase, which 
will complicate the mapping process. Note, further- 
more, that only the planet's permanent day-side could 
be mapped using reflected light, in such a case. To con- 
strain the latitude-dependence of albedo, one would need 
a high-obliquity planet and observations spanning many 
different phases. For our analysis we use two classes 
of model maps: one constructed from sinusoidal vari- 
ations in albedo as a function of longitude, and the 
other with uniform l ongitudinal slices of constant albedo 
() Cowan fc Agoll [2008 '1 . In both cases the albedo is con- 
stant with latitude. We compare these two models ex- 
plicitly in Appendix II. 

The transformation from to A(0) to A*{t) is essentially 
a low-pass filter which preferentially preserves informa- 
tion about large-scale color variations on the planet. The 
smoothing kernel (the product of V and /) has a FWHM 



6 



Cowan et al. 



of 76° for the March observation and 67° for the June ob- 
servatfon. Furthermore, the transformation can be par- 
ticularly insensitive to certain modes and these cannot 
be recovered from the light curve inversion (for example, 
a planet obse rved at full phase has invisible od d sinu- 
soidal modes, |Ruie3[l90& £(ManXAgoi [200^ Even 
in the idealized toy example discussed in Appendix II, 
therefore, the deconvolution can only recover the broad- 
est trends. 

These best-fit parameters and their uncertainties 
are determined using a Markov Chain Monte Carlo 
(MCMC), opti mized to h ave acceptance probabilities 
near 25% (e.g.. lFordll2005f ). The uncertainty in the si- 
nusoidal model coefficients tends to increase with the 
harmonic index n because high-frequency modes have 
a relatively small impact on the observed light curves. 
We truncate the series when we have sufficient terms to 
get a reduced of order unity. In this study we find 
that including modes up to rt = 3 or 4 is sufficient, de- 
pending on the phase angle and the color being used. 
By the same token, 7-slice and 9-slice models are used 
in our fits since they have the same number of free pa- 
rameters. The sinusoidal and N-slice longitudinal maps 
of the red and blue eigencolors (roughly, land coverage 
and ocean coverage) we construct using PCA and light 
curve inversion are shown in Figures [6H9l The March 
and June maps differ due to different cloud cover on the 
two days. Nevertheless, the broad peaks and troughs 
occur at the same longitudes at both epochs, indicating 
that the eigencolors are sensitive to permanent surface 
features, not merely clouds. The red eigencolor is more 
sensitive than the blue eigencolor to the positions of con- 
tinents and oceans on Earth, despite the fact that both 
eigencolors can be fooled by clouds. This is because the 
red eigencolor is most sensitive to the near-infrared light 
that arid regions of Earth reflect at, and those regions 
are generally cloud-free. 




longitude [degrees] 

Fig. 6. — Comparison of N-slicc (blue) and sinusoidal (black) lon- 
gitudinal maps of the red eigencolor based on the March EPOXI ob- 
servations (top panel) to the MODIS cloud-free, equator-weighted 
distribution of continents (bottom panel). The thickness of the 
lines in the top panel show the la uncertainty on the maps. The 
red eigencolor is particularly sensitive to arid regions (the South- 
West USA, the Sahara and Middle East, and Australia) since they 
are consistently cloud-free. 

For our baseline model we assume that the planet's 
rotation axis is perpendicular to its orbital plane (zero 



0.06 




longitude [degrees] 

Fig. 7. — Comparison of N-slice (blue) and sinusoidal (black) lon- 
gitudinal maps of the red eigencolor based on the June EPOXI ob- 
servations (top panel) to the MODIS cloud-free, equator-weighted 
distribution of continents (bottom panel). The thickness of the 
lines in the top panel show the Itr uncertainty on the maps. The 
red eigencolor is particularly sensitive to arid regions (the South- 
West USA, the Sahara and Middle East, and Australia) since they 
are consistently cloud-free. 




longitude [degrees] 

Fig. 8. — Comparison of N-slice (blue) and sinusoidal (black) 
longitudinal maps of the blue eigencolor based on the March 
EPOXI observations (top panel) to the MODIS cloud-free, equator- 
weighted distribution of oceans (bottom panel). The thickness of 
the lines in the top panel show the Icr uncertainty on the maps. The 
fit to Earth's latitudinaly averaged water content is not excellent, 
notably the Atlantic Ocean and the Americas. The discrepancy is 
due to cloud cover, as discussed in the text. 

obliquity) and determine the best-fit longitudinal map 
in each of the two principal colors, assuming the same 
underlying map for the March and June observations. 
The June map of the red eigencolor, the most impor- 
tant of the principal components, is shown in Figure [10] 
compared to a cloud- free MODIS^ map of land coverage. 
Cloud cover, shown in Figure [TTJ keeps the match from 
being perfect, but our blind analysis of the light curves 
clearly picks out the Atlantic and Pacific Oceans, as well 
as the major landforms: the Americas, Africa, and Asia. 

4.3. Obliquity 

The obliquity of exoplanets, the angle between their ro- 
tational and orbital axes, will not be known a priori. The 
simplest assumption — and the one used in our baseline 

http://modis.gsfc.nasa.gov/ 
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longitude [degrees] 



Fig. 9. — Comparison of N-slice (blue) and sinusoidal (black) lon- 
gitudinal maps of the blue eigencolor based on the June EPOXI ob- 
servations (top panel) to the MODIS cloud-free, equator-weighted 
distribution of oceans (bottom panel). The thickness of the lines 
in the top panel show the la uncertainty on the maps. The fit 
to Earth's latitudinaly averaged water content is not excellent, no- 
tably the Atlantic Ocean and the Americas. The discrepancy is 
due to cloud cover, as discussed in the text. 




Fig. 10. — AitofT projection showing the land distribution on 
Earth in a cloud-free MODIS map (top panel) and the distribution 
of land as determined from the June disc-integrated EPOXI light 
curves (bottom panel). The EPOXI map has a longitudinal reso- 
lution of approximately 60°; it has no latitudinal resolution, but is 
weighted toward the equator due to viewing geometry. 

model — is zero-obliquity, but this cannot be assumed 
as a rule. We fit the light curves using different assumed 
orbital axes and verify that the longitudinal maps are in- 
sensitive to this assumption. For example, the longitudi- 
nal maps obtained assuming the Earth's correct obliquity 
of 23.5° are almost indistinguishable from our baseline 
(no obliquity) case. This is not surprising. The appar- 
ent albedo of the planet amounts to a weighted average 
albedo of the regions of the disc that are both visible and 
illuminated. During a single rotational period all longi- 
tudes of the planet will be both visible and illuminated 
for some fraction of the time. The albedo variations on 
an oblique planet must be larger to account for the same 
observed light curves, however. We see this effect in our 
models: the longitudinal maps for the 23.5° obliquity 
have slightly greater amplitude variations than those for 
the zero obliquity case, while in the case of 90° obliq- 
uity the subtle changes in apparent albedo in Figure [1] 



Fig. 11. — Aitoff projection showing the fractional cloud coverage 
for 03/18/2008 (top panel) and 06/04/2008 (bottom panel), based 
on MODIS Aqua. The regularly spaced white artifacts represent 
missing data. 



can only be explained by enormous (100%) changes in 
albedo from one region of the planet to another. 

Obliquity also determines which regions of the planet 
can be mapped. As long as the planet's rotation axis 
resides in the sky plane, the reflected light we observe 
is preferentially coming from near the equator. If in- 
stead the planet's rotational axis is not in the sky plane, 
we preferentially "see" some non-zero latitude. A lon- 
gitudinal map can only be a good representation of the 
latitudes that are both visible and illuminated, and the 
latter will change depending on the phase of the planet. 
During the March EPOXI observations, the sub-solar 
and sub-observer points were both close to the equa- 
tor; for the June observations the sub-observer was again 
nearly equatorial but the sub-solar point was at 22° N 
of the equator (it being northern summer) . The minus- 
cule effect of this change in viewing geometry can be 
seen in Figure [121 The differences between the March 
and June observations must therefore be due to different 
cloud cover. 




180 w 



90 w 

longitude [degrees] 



90 E 



180 E 



Fig. 12. — The effect of a change in viewing geometry: the 
equator-weighted longitudinal water distribution on Earth based 
on the MODIS map is shown in the solid line. The dotted line 
shows the same but weighted in favor of 11° N, appropriate for 
a viewer above the equator near summer solstice. The effect is 
effectively negligible, justifying our assumption of zero obliquity. 
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5. DISCUSSION 

Spectra of habitable terrestrial exoplanets will tell us 
the structure and composition of their atmosphere, but 
may require integration times of weeks to months. Since 
this is longer than the rotation period of most plan- 
ets, spectroscopy can only tell us about the spatially- 
averaged planet. Photometric light curves, with integra- 
tion times of hours to days, have the potential to reveal 
sp atial variations in the planet's properties. Simulations 
bv lFord et al.l ()2001[ ) indicated that diurnal variability in 
the albedo of an unresolved, cloudless exoplanet could be 
used to determine its ocean versus land fraction. But in 
their models it was the specular nature of oceans, rather 
than their blue color, that distinguis hed them from con- 
tinent s. The more detailed work of IWilliams fc Gaidod 
(|2008l ) indicates that on a cloudy world like Earth, the 
contribution of specular reflection from oceans will be 
tiny compared to the diffuse reflection from clouds. The 
models of Palle et al.l HQOS) showed that despite changes 
in cloud cover, the diurnal albedo variations of an Earth- 
like planet could be used to determine it's rotation rate. 
Our study shows that with o bservations q ualita tively 
similar to those considered by iPalle et "all ()2008D . but 
with greater signal-to-noise and better spectral resolu- 
tion, it is possible to actually map the longitudinal distri- 
bution of colors — and by extension the dominant surface 
types — of Earth. 

Earth clouds have higher albedo at all seven wavebands 
than ocean or continents (see first the bottom panel of 
Figure [3]) , a nd most regions o f the planet have variable 
cloud cover (jPalle et al.ll2008l ). Observations of multi- 
ple consecutive planetary rotations would yield multiple 
similar longitudinal maps. If differences between these 
maps were attributed to changes in cloud cover it would 
be possible to create maps of average cloud cover at each 
longitude, as well as maps of cloud variability. It may 
even be possible to partially "remove" clouds since the 
lowest albedo at each longitude would correspond to the 
observation with the least cloud cover. It would be im- 
possible, however, to strip the clouds from regions that 
are permanently shrouded (e.g.: tropical rain forests). 
Insofar as such clouds are permanent features, they can 
be mapped as terrain features, like oceans and conti- 
nents. 

Clouds are doubly important because they change over 
time and dominate the total albedo of Earth. However, 
insofar as roughly half of Earth is cloud-covered at any 
point in time, the bolometric albedo of the planet does 
not change much over 24 hrs (as shown in bottom right 
panel of Figure 1). What can change more significantly is 
the reflected color of the planet, which is what we have 
studied in this paper. Since clouds are roughly gray, 
any color (apart from blue Rayleigh scattering) is com- 
ing from cloud-free regions. What we have shown in this 
paper — without prior assumptions — is that these cloud- 
free regions come in two varieties: blue and red. It is our 
contention that these correspond to oceans and conti- 
nents. Moments when both the red and blue eigencolors 
are high correspond to moments when a particularly high 
fraction of the visible, illuminated Earth was cloudy. It 
is incorrect to think that the color variability is simply 
a measure of cloudiness, however. Had that been the 
case, the PCA would only have found a single significant 



eigencolor, whereas it found two. Despite the fact that 
clouds dominate the over-all albedo of Earth, the pres- 
ence of continents and oceans is still discernible in the 
disc- integrated colors of the planet. 

Note that a blue broadband spectrum does not — in 
and of itself — imply liquid water on the surface of a 
planet (e.g., Neptune). The spatial variability in the blue 
color is significant, however. An alternative explanation 
for spatially inhomogencous blue colors could be partial 
cloud coverage: the increased path length in regions with 
fewer clouds could increase the importance of Rayleigh 
scattering. But the spotty cloud cover on such planet 
might, like time- variable cloud cover, give away the pres- 
ence of water near its condensation point. Furthermore, 
the blue patches on such a planet might reveal themselves 
by their very steep blue spectrum. In conjunction with 
time-averaged spectra, broadband light curves therefore 
provide a powerful test for the presence of liquid water 
on a terrestrial planet, and hence habitability. 

6. CONCLUSIONS 

In this paper we have shown that despite simplify- 
ing assumptions (edge-on geometry, zero-obliquity, dif- 
fuse reflection) and the presence of obscuring clouds, one 
can use time-resolved photometry to detect and map vast 
blue surfaces, separated by large red regions. If we saw 
such features on an extrasolar terrestrial planet it would 
strongly suggest the presence of continents and oceans, 
and indicate that the planet was a high priority for spec- 
troscopic follow-up. 

Although Earth is most reflective at short wavelengths 
due to Rayleigh scattering, the wavelengths longward of 
700 nm provide the most spatial information because the 
relative diurnal variability at those wavelengths is great- 
est (25 -30%, in qualitat ive agreement with the simula- 
tions of lFord et al.ll200ll ). as shown by the error bars in 
Figure [3l The relative variability is greater for smaller 
illuminated fraction, in accordance with geometric con- 
siderations. The logical conclusion is that observations of 
exoplanets at crescent phase would provide the greatest 
spatial resolving power, but this is not true in practice 
due to poorer signal-to-noise ratio: the flux from a cres- 
cent planet is smaller than at quadrature and at small 
angular separations the planet is lost in the glare of its 
host star. Note that interesting measurements can be 
made if the planet passes directly in front of (transit) 
or behind (secondary eclipse) its host star, but for habit- 
able terrestrial planets the odds of this are not very good 
(e.g., an Earth-analog has a 0.5% probability of transit- 
ing a Sun-like host star). 

The unknown obliquity of exoplanets does not rep- 
resent a serious obstacle to mapping their longitudinal 
color variations, although it will affect which parts of 
the planet are being mapped. There are pathological or- 
bital configurations that would prohibit mapping (e.g., 
pole-on rotation axis observed at full phase) but planets 
in such configurations will be impossible to observe in 
any case. Planets will typically be observed near quadra- 
ture, so longitudinal maps will not depend sensitively on 
obliquity. 
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APPENDIX I: VISIBILITY & ILLUMINATION 

One can define a right-handed orthonormal coordinate system, x, y, z, in the observer's inertial reference frame with 
the origin at the center of the planet, the a;-axis extending towards the observer, and the y- and 2;-axes in the sky 
plane. The orbital and rotational angular velocity vectors of the planet in this frame are Worb and Wrot, respectively. 
If the planet is in an eccentric orbit, the amplitude of oJorb will be a function of time, but in any case the direction 
of the vector is constant. In the interest of simplicity we consider a circular orbit. Neglecting precession, w^t will 
be a constant vector (note that this is the rotation of the planet in an inertial frame, not the rotation with respect 
to the host star). We define the position of the sub-stellar point at f = as ^star and the intersection of the prime 
meridian and the equator on the planet as fequ- Note that fstar is always perpendicular to Uoib, and iequ is always 
perpendicular to Wrot- 

At t = we define the orthonormal coordinate system fixed with respect to the planet surface uq = iequ) wq = Ld^ot 
and Wo = m)o X uq. Likewise, we define at t = the orthonormal vectors defining the star reference frame: do = ^star, 
fo = Worb and eo = /o x da. 

The position of a region on the planet can be described by its complementary latitude, 9, measured from the planet's 
north pole, and its east longitude, (f), measured along the equator from the prime meridian. This patch has a position 
in the observer frame f = sin cos (/'■u -|- sin sin t/)?) + cosOw. The visibility of this region from the perspective of the 
observer is V{6, (j), t) = max[f • 0]. The inner product can be expanded as: 

f ■ X = sin ^ cos((/> -|- LVrott)uo • x 

+ sin 6 sm{(j) + u>rott)vo ■ x (4) 

-|- COS^Wo • X. 
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The illumination of this patch is I (9, (f), t) — max[f • d, 0]. The inner product can be expanded as: 

f ■ d = sin 9 cos{LUoiht) coa{if) + WrotO'^^o • do 
+ sin 9 cos(ijJorbi) sin(0 + cjroti)wo • do 
+ sin 9 sin(worbi) cos(0 + uJiott jiio ■ eg 
+ sin 9 sin(a;orbi) sin((/) + LJrott)vo • eo 
+ cos0cos(a;orbt)wo • do 
+ cos0sin(a;oibi)w;o • eo- 

The visibihty and illumination can be expressed compactly in terms of the sub-observer longitude, 0obs(O — '/'obs(O) — 
uJrott, and the constant sub-observer latitude, ^obs, as well as the sub-stellar longitude, ipstar, and latitude, flstai- 

V — max[sin 9 sin ^obs cos(0 — (pohs) + cos 9 cos ^obs , 0] 

1= max[sin6'sin6'star C0S((/) - (/)star) + cosflcos^starjO], 

where the sub-stellar longitude is related to the orbital phase measured from the solstice, ^{t) = + Worbi, and the 
constant planetary obliquity, 6'obi, by cos6'star = cos^ sin6'obi- 

APPENDIX II: COMPARING N-SLICE AND SINUSOIDAL MAPS 

The N -Slice and sinusoidal models are compared in Figure 1131 with the resulting light curves shown in Figure [T3] 
(see also lCowan &: Agoill2008D . A MODIS map of liquid water content (s.f. top panel of Figure [TU|) was integrated to 
a one-dimensional, equator- weighted map of water content, the black line in Figure [TSl A model light curve, shown in 
black on Figure [TH was generated assuming photometric uncertainties of 1%. This light curve became the input for 
light curve inversions using sinusoidal (red) and N-slice (blue) maps, shown here with ±1(7 intervals. 



Actual Map 




longitude [degrees] 

Fig. 13. — Comparison of sinusoidal and N-slice longitudinal maps. The thickness of the lines denotes the ztlcr intervals. The basic 
features of the underlying map (the major continents and oceans) are recovered by either mapping technique. 
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Fig. 14. — Comparison of the light curves produced by the three maps of Figure [131 
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APPENDIX III: PLANET MAPPING IN THE SPECULAR REGIME 

An interesting limiting case arises when the entirety of the light from the planet originates fr'om the glint spot where 
the product VI is maximized, as would be the case for purely specular reflection. The latitude of the specular point 
is given by 

cos 6>star + COS 6>obs 
^2[1 + Sin6lstar sin 6lobs COs((/)star - 0obs) + COS 6'star cos 6lobs] 

and its longitude is given by 

sin 0star Sm «star + Sm 0obs Sm 6>Qbs 

cos 4>ohs sin 0ohs 

In Figure [T5] we show the results of specular inversion on the light curves shown in Figure [S] We have assumed 
edge-on, zero-obliquity geometry. The map of the red eigencolor shows three distinct peaks, corresponding to (from left 
to right) North America, Africa and Asia. The map of the blue eigencolor is characterized by a high plateau between 
90° E to 90° W, which corresponds to the Pacific Ocean. The diffusely and specularly reflecting cases bracket more 
realistic scattering phase functions, with Solar System planets and moons being closer to the former. The detectability 
of the major continents and oceans on Earth using either assumption indicates the robustness of the result. 



cos tlspoc = . . ^ . ^ . , , X . ^ (^) 



Dill v^star c'iii i^star ^ ^^^'^ •Tabs ^^'^^ ^obs /c\ 

tan0spoc = 7 r-7 , ^ ZHT^ — ■ i^) 




Fig. 15. — Longitudinal maps of land (red eigencolor) and water (blue eigencolor) on Earth, based on a principal component analysis of 
disc-integrated light curves and the assumption of specular reflection. 



